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FINITENESS OF THE COULOMB GAUGE QCD 
PERTURBATIVE EFFECTIVE ACTION 

A. ANDRASI AND J. C. TAYLOR 


Abstract. At 2-loop order in the Coulomb gauge, individual Feynman graphs 
contributing to the effective action have energy divergences. It is proved that 
these cancel in suitable combinations of graphs. This has previously been 
shown only for transverse external fields. The calculation results in a general¬ 
ization of the Christ-Lee term which was inserted into the Hamiltonian. 
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1. Introduction 

In QCD, the Coulomb gauge has some attractions. It is the only manifestly 
unitary gauge (no ghosts). It may be convenient in practice when there is a natural 
rest frame, as in thermal QCD or for heavy quarks. It has been used in studies of 
confinement (see for example 0). It has been used in lattice calculation, see for 
example m), and in studies of infra-red behaviour (see 0)- But, in perturbation 
theory, it has complications due to energy divergences, that is Feynman integrals 
which are divergent over the internal energy integrals (for fixed spatial momenta). 
These divergences appear at 1-loop order, but in this case are easily curable (see 
most simply by using, instead of the Lagrangian formalism, the phase-space, 
Hamiltonian formalism, as we do in this paper. 

At 2-loop order, there are more subtle energy divergences. These occur in diver¬ 
gent integrals of the form 
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(capital letters, P,Q,R etc denote spatial momenta and P = |P| etc). This prob¬ 
lem is resolved if it can be shown that graphs can be grouped so that integrals like 
(1) appear in the combination 
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( 2 ) 


= -(7rV3)F(P,Q,R) 


as may be shown by contour integration. 

For definiteness, suppose that we are calculating, to 2-loop order, the effective 
action r(A“) as a functional of the external vector field A“ (i, j,.. are spatial indices 
and a, h,c,.. denote colour). (Parts of F depending on Ag or Ef contain no divergent 
integrals like (1).) 

The requisite theorem, that graphs can be combined to give (2), was proved by 
Doust [B] (see also [7]) for the special case of transverse external fields, dtA^ = 0. 
But the effective acion is formally defined (see the Appendix) for general A“, and 
this is indeed required in order to state the BRST identities for F. The purpose of 
the present paper is to extend Doust’s theorem to general A“. The 0{g^) case was 
already studied in [5]. The divergent integrals (1) are connected with the fact that 
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operator ordering of the Hamiltonian (which is non-polynomial) may not be well- 
defined. Christ and Lee [5] (and references therein) argued that the correct ordering 
entailed adding onto the naive Hamiltonian two operator functionals, called Vi and 
V 2 , of the operator which in the Coulomb gauge is transverse by definition. To 
avoid double counting, the Feynman integrals (1) and (2) would have to be omitted. 
Doust found that the integrals (2) contributed the same functionals Hi -|- V 2 to the 
effective action. 

1.1. The flow gauge. In order to manipulate integrals like (1) rigorously, a tem¬ 
porary regularization is necessary. Dimensional regularization (which we use, with 
space-time dimension d, to control ordinary ultra-violet divergences) does not have 
any effect on (1). So we (following Doust) employ a ’flow gauge’ with a parameter 
0. For finite 0 the function F in (1) depends also on 0po,0qo,0ro so as to make 
the integral converge. So in (1) and (2) we must replace F by (we write p for the 
Lorentz vector (po; P)) 

(3) F(P, Q, R, 0po, 0qo, 0ro) = F(p, q, r). 

The Coulomb gauge is the limit 0^0, which we take only after getting the 
combinations (2). (0 = 1 is the Feynman gauge, so the flow gauge is sometimes 
called an interpolating gauge). The flow gauge is defined by the gauge-fixing term 

(4) 

(the limit 0 —5> 0 imposes the transversality condition on the operator H“). The 
effect of this is that the bare Coulomb propagator, — 1/P^, is everywhere replaced 
by 

(5) 1 ^ 

^ ^ p2 _ 02p2 - -p2 

proof.tex This is what causes the function F in (1) and (2) to be replaced by (3). 

In the Coulomb gauge, ghost loops exactly cancel closed Coulomb loops. But in 
the flow gauge this is not exactly the case. Ghosts, but not Coulomb loops, have 
a coupling O{0^) to the Coulomb field. Hence, as well as graphs similar to Fig.l, 
graphs with ghost loops have to be included in the proof. When the convergent 
form (2) is finally attained, the ghost graphs can be omitted. 

The method of our proof, following [6], is to show that F{p,q,r) is a cyclically 
symmetric function of its arguments. Then the integrand in (1) can be symmetrized, 
to give (2). Finally, having attained the convergent combination in the square 
bracket in (2), the limit 0—^0 can be taken in F in (2) giving the Coulomb gauge. 

1.2. An example. An example (at 0{g^)) of a Feynman graph leading to an 
integral of the form of (1) is shown in Fig.l. Our notation for Feynman graphs is 
explained in Fig.2. Dotted lines denote the propagator (5). Dashed lines denote the 
propagator for A“ (transverse in the Coulomb gauge, but not transverse for nonzero 
0). Solid lines always denote a momentum (either spatial or time-like), multiplied 
by i, in the numerator (the sense of the momentum being directed towards the 
neighbouring vertex). The external dashed lines represent the external field A“. In 
Fig.l, the two denominators in (1) come from the lines labeled p' and q in Fig.l. 

In the Coulomb gauge {0 = 0) all the relevant graphs have a ‘spine’ coming 
from the Coulomb operator in the Hamiltonian. In Fig.l this spine consists of the 
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Fig.1 An example of a graph O (g^) giving an integral of the form (1). 


p,r,ri,r',q',qi lines. In addition to this spine, there are two lines {q,p') where 
Ai and Ej operators in the Coulomb operator have been contracted with one- 
another. In the flow gauge, the structure of the graphs are the same though their 
interpretations are not as simple. 

In Fig.I there is just one line, labeled ri, which has a dotted line with no solid 
line on it. In all the relevant graphs there is on the spine just one line like this, 
which may be in any position, on the centre line or the righthand or lefthand one. 
(To order (/” there are n — I such possible positions). In the case of transverse 
A“, treated by Doust, the contribution is independent of the position of this line, 
which simplifies the argument. We show that the proof can be made without this 
simplifying property. 

We need a compact notation in order to write down the general graphs giving 
integrals of the form of (1). This we explain in the next section. 


2. The ghost and Coulomb propagators 

The example in Fig.l has three chains of lines. Each chain consists partly of a 
propagator in the external field Af. On the left and right in Fig.l, these propagators 
have the same form as the ghost propagator (though they are not ghost lines), and 
in the centre the Coulomb propagator occurs. We need a compact notation for 
these two propagators, which we now give. 





QCD EFFECTIVE ACTION 


5 


j 


Err. 


m 

n 

Em 

Aq 

0 

n 


(-1) 

(27r)^i 

1 

{2'n-yi 


1 


+ it] 




KrK, 


- 0‘^kl 

1 


K'^ - 9'^ kl 


iko 

WTiv r™^ ~ K‘^-e‘^kl\ 




k'^ + irj 


EmEn 




iErr 


K"^ - e'^kl 

iErr 

^ K'^-9‘^kl 
1 

^ K'^-e'^kl 


for each closed ghost loop 
for each vertex 
for each propagator 



. b, 0 

c,k ■ a,m 

- gf-»Cg^, 


Fig. 2: The graphical notation for Feynman propagators and the rules for ver¬ 
tices in the flow gauge. Doubled lines represent ghosts. 


We call the ghost propagator G°'°' (P, P'; Opo, Op'^). In the Coulomb gauge, there 
would be no dependence on porP'o, but in the flow gauge, such dependence is intro¬ 
duced through the denominators (5). 
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G“{p.fO 
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C*(p.p-) 

I I 




... 


Flg.3 Definitions of graphical symbols and relations between them. 
The first two definitions oonsspond to equations (7) and (9). 


Define the antisymmetric matrix field 

( 6 ) Af{k) = r^^At{k) 

where A^(k) = A'^(—k) is the external field in momentum space. 

The ghost propagator G is defined by the integral equation 

(7) -p2G‘^‘^'ip,p') = S--'S‘^ip-p') + zgI 

This may be solved iteratively to give a series vci g. As an example, the second 
order term is 

(8) g^ J d'^pi\P^lfP^]-^V.A^\p^ -p)Pi.A^“'(pi -p'). 

The Coulomb propagator (7““ (jp,p') is given in terms of the ghost propagator 

by 

(9) G““'(p,p') = J d’^PiG^\p,Pi){-^l)G^'\-p'-pi). 

As an example, one of the three second order terms in G is 

g^ J d''Pi[p2p2p72]-ip -p)P'.A“'*'(p' -pi). 


( 10 ) 
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The example (10) occurs on the centre chain in Fig.l (with p replaced by r). 
Note the difference between (8) and (10) in general, but for transverse A they 
are equal. 

The series for G and C are illustrated graphically in Fig.3. 



Fig.4 The structjre of graphs contributing to the divergent integral (1). 
Each line represents a chain cf arbitrarily many propagators, 
as defined in Fig.3. 


3. The energy-divergent graphs 


We need a notation for graphs in which an arbitrary number of external gluon 
fields are attached to the three chains. This notation is shown in Fig.3. Then the 
six types of graph which give an integral containing (1) (with the substitution (5)) 
are shown in Fig.4. They are all drawn with the p-chain on the left and the g-chain 
on the right. These graphs give contribntions of the form 
( 11 ) 


^( 277 )- 


d jabc jao c 


d‘^pd\d\S{p + q + r)d‘^p'dU'd‘^r'6{p' + q' + r') 


p2 q/2 


B, 


where the factors B will be expressed graphically in Figs.7 and 8 later in this 
section. (The factor 1/2 in (11) is needed because the set of graphs which we use 
has symmetry under the interchange of the two internal vertices.) 



8 


ANDRASI AND TAYLOR 


In order to connect (11) to (1), we need a lemma. This states that, up to terms 
0{9), we may replace or p'^qo/(p'"^q"^) by 

Po Qo 


( 12 ) 


p^ 


le 


9o 


■ le 


(reminding ourselves of the ie factors, which we usually leave understood). The 

proof is as follows. The difference is an integral of the form 

(13) 


dpodqodro6{po + qo + ro) 


PoQo Poqo 


p2q/2 


p2q2 


F(P, Q, R; Opo, Oqo, K^™)) 


where the {m = 1, 2, 3..) are external gluon momenta. Also q^ — q^ is equal to 
a linear combination of a subset of the . We now make a change of variables 
(14) Po = OpQ, qo = 9qo,ro = 9ro. 

Then, in terms of the new variables of integration, the difference between the two 
terms in the square bracket in (13) lies in dependence on 9^P'^,9^P’^,9k^\ The 
integral is convergent, because of the terms in F like 

1 


(15) 


- 6»2R' 


2 ’ 


so we may let d —>■ 0, and the integrand in (13) is zero in the limit. 

In view of this lemma, the square bracket in (11) may be replaced by (1). Having 
done this, the contributions to B in (11) may, with the use of the definitions in 
Figs.3 and 5, be written as in Fig.7. All these graphs have at least one external 
field attached to the p-chain and at least one to the q-chain (these are the external 
dashed lines shown explicitly; in addition there are arbitrarily many attachments 
in the C and D propagators). But there are the special cases in which there is no 
external attachment to the p-chain, or none to the q-chain, or none to either. Some 
examples of these cases are shown in Fig.8. Our aim is to show that the sum of 
the graphs in Figs.7 and 8 has symmetry under cyclic permutations of 

(16) (p,p';a,a'), {q,q]h,b'), {r,r,c,c). 

Then, with the integration (11), (1) can be replaced by the convergent integral (2). 

To order p", there are contributions with either n-|-lornorn — 1 denominators 
like (5). It turns out that each of these three classes separately has cyclic symmetry; 
so we treat them one by one in the following subsections. 


3.1. (n -I- 1) denominators. For these terms, in the graphs in Fig.7, we take the 
PiPj/P"^ and p 2 /p 2 terms from Fig.5. We then use the first two identities in 
Fig.6, which follow from (7) and (9), thus getting one contribution from (a) and 
two each from (b), (c), (d), (e) and (f). In the ghost graphs (d), (e) and (f), we 
use also the third identity in Fig.6, which follows from momentum conservation. 
As a result of these steps, there are finally four contributions from each of (d), (e) 
and (f). The total number of contributions from all six graphs is thus 17. There 
are sign factors. The ghost graphs (d), (e), (f) each have a minus sign. There are 
additional minus signs for (a), (d), (e) and (f) because there colour matrices are in a 
different order from that in equation (11). Finally, contributions from the last term 
in the second identity in Fig.6 receive another minus sign. The term from (a) is 
canceled by one of the contributions from (d), so that fifteen graphs remain. These 
fall into five triplets, each triplet being invariant under cyclic permutations of the 
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Fig .5: Definitions of graphical symbois and relations between them. 


three chains, that is cyclic permutations of the variables in (16). This allows the 
divergent integral (1) to be replaced by the convergent integral (2). Three of the 
triplets are shown in Fig.9. The remaining two are obtained simply by reflection, 
that is the change of variables 


(17) 


P,P, q, q, r, r'] a, a', b, b', c, c fA q, q ,p,p , r'r; b, 6', a, a, c, c. 


In some of the graphs in Fig.9 we have made use of the change of variables 
(18) p, q, r; a,b,c<r^ -p', -q', -r'; a', b', c' 


that is inversion of the graph. When this has been done it is indicated by a prime, 
e.g. (c) ^ (c)' 

In special cases, some graphs from Fig.8 have to be added to those from Fig.7, 
in order to contribute the last term in the first identity in Fig.6. 


3.2. n denominators. For these terms, we take from Fig.7 just one term with 
Sij or with I. Because of symmetry under (18) (inversion), we can wihtout loss of 
generality take these to be at the top of the graphs. Then the possible graphs are 
shown in Fig.10, with the notation in Figs.3 and 5. At the top vertex, momentum 
conservation takes the form 


( 19 ) 


d‘^kS'^{p + q + r-k)Af{k) 
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Fig.6 Three Identities expressed graphlcally.The first two follow from 
the def ning equations (7) and (9) or from Flg.3. The third is a 
consequence of momerrtum conservation P Q R = 0. 


where (k) is the external field shown explicitly by the dashed lines in Fig. 10. 
Expression (19) is symmetric in p,q,r. 

In Fig. 10, use has already been made of the first identity in Fig.6. To simplify 
these graphs, we use the third identity in Fig.6 for graphs (d) and (f), and apply 
Jacobi identity, shown graphically in Fig.11, to graphs (d), (ei) (eii). This identity, 
stated algebraically, is 

( 20 ) 

where is defined in (6). There result 14 graphs. Graph (a) cancels one of 
the contributions from (d), when we make use of (17). Again using (17) where 
necessary, we obtain 12 graphs, which fall into four groups of three, each group 
having cyclic symmetry. These are shown in Fig.12, labeled with the names of the 
graphs in Fig.10 from which they originate. Again, the cyclic symmetry enables us 
to replace the divergent (1) by the convergent (2). 

3.3. (n — 1) denominators. In this case, we take two terms Sij or 1 from Fig.5 in 
Fig.7. The results are shown in Fig. 13. Then we use the Jacobi identity Fig. 11 twice 
in graph (d) and once each in graphs (e) and (f). There result 11 graphs. Graph 
(a) cancels one of the contributions from (d), leaving nine graphs. These fall into 
three cyclically symmetric triplets, as shown in Fig. 14, having used the inversion 
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Fig.7 Graphs daivad from Fig.4. which contributs to the factor B 
In equation (11). The notation Is defined In Figs. 3 and 5. 


symmetry in (18) where necessary. Once again, the cyclic symmetry enables the 
divergent (1) to be replaced by the convergent (2) in (11). 

This completes the proof for graphs in which at least one external field A is 
attached to each chain. But there are exceptional graphs for which this is not the 
case, like those in Fig.8. We must now deal with these exceptional cases. 

4. Exceptional graphs 

These are graphs with either no external gluon attached to the g-chain or none 
attached to the p-chain. Because of reflection symmetry (17), we can without loss 
of generality take the former case. The graphs are shown in Fig.8. Again we have 
to consider the cases with (n + 1), n, or (n — 1) denominators, but the first of these 
is already included in the previous section, so we need only consider the latter two 
cases. 

4.1. n denominators. The relevant graphs are shown in Fig. 15, where the central 
rectangle denotes S’’’’ in the graphs (d), (e), (f); and S’’’’ Sij in (a), (b), (c). The 
graphs (cf), (cii), (ei), (eii) arise because of the use of the second identity in Fig.6. 
The last identity in Fig. 6 (momentum conservation) is applied to graph (d) in Fig. 15 
to give the three graphs in Fig. 16. 
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Fig.B Grap^ 1 s derived from Fig.4. which contribute to trie facrior B in 
equation (11), and which have no external field attached to the 
q-line. The labels correspond to those In Fig. 7. 


We then prove symmetry under the interchange of the p- and r-chains. Graph 
(a) in Fig.15 cancels graph (dii) in Fig.16. Graphs (cii) and (eii) in Fig.16 each 
separately have the required symmetry, making use of inversion invariance under 
(18). And the two pairs 

(21) {(ci)', (dy)}; {{ei),{da)} 

each are invariant under the interchange. 

This means that in (1), 

(22) F{p,q,r) = F{r,q,p) 

(i^(P,Q,R) becomes F{p,q,r) in the flow gauge, as in (3)). Therefore we can 
replace (1) by 

(1/2) [ dpodqodro6{po + qo + ro)F{p,q,r) 


(23) 


Po 


do 


do 


Pq- do ~ do - + ie b 


We note also that, in Figs.(15) and (16), F factorizes 
(24) F = <l>f iF,0VoW {^,0^1), 


A)_ 

+ ie 
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Flg.9 TOroo tlplets at s^p^8 obtained from Fig . 7 by use of tbe Mentiles In Figs. 3 and 
6. Eadi triplot is invariant mder cyclic pemulaticrs of its tlwee ctiahs. In addition 
Ihere are tvro mne biplalB oblaned from the first and second tiipteta above by the 
liansfciiTiatlon (17) (lef acflon) Tin IstteiB under the graphs Indicate the ^aptis tii 
Flg.7 frori vrhleb they ere derived. TTn ptiine on (e)f etc. Indicates that the 
barsfarmation (16) (hverston) has been made, lha aigns attached to these gr^ihs 
come about as explaired in ttie tout 


neglecting terms 0(9k\^^) where are external momenta; so 

(25) J = 0{e), 

and similarly for the r-integral (just as at one loop order). This allows us to insert 
the third term (poTo/p^t^) into the square bracket in (23) and then use the identity 
in (2) to get 

(26) - {^^/2)F 

and finally let 0 —>■ 0. Thus we have achieved our aim, but there is the important 
difference from section 3 that we here have a factor of a half rather than a third as 
in (2). 

4.2. (n — 1) denominators. The relevant graphs are in Fig. 17. For graph (d), we 
use the Jacobi identity (Fig.11) twice to get the three graphs in Fig.18. As usual, 
graph (a) is canceled by (dl)’. The remaining four graphs fall into two pairs 

(27) {(e), (d3)}; {(c)', (d2)} 

each of which has symmetry under the interchange of the p- and r-lines. After that 
the argument continues exactly as in subsection 4.1. 
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Fig. 10 Grep^s obtained from Flg.7 vritb n denomktaton. The rotaton Is explalred In 
Fig.5. The second identity in Fig.ti has been used, giving (d), (csjend (ei)i (di). 
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Fig. 11 The ^aplical Jacobi identity expresdng equator ( 20 ) in tire text. 


There is a special case in which, in Fig.4, both the p-chain and the g-chain in 
graph (a) consist of the single terms dij, and in graph (d) of the single factors 1. 
The result is just the same as in [^, and is 

(28) Tr^{d-2) + Q + 
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which they are dailved The double prime denotes the transtbrm by (17) (iat1ectlan)i. 


In dimensional regularization this integral is zero, but we write it here for complete¬ 
ness. Below we neglect all such terms which vanish in dimensional regularization. 


5. The generalized Christ-Lee functions 

In both section 3 and section 4, the energy integral has been shown to be finite, 
and could be done using equation (2). But in section 3 this resulted in a factor 
—7r^/3 while in section 4 the factor was —7r^/2 (leaving aside (28)). The limit 6—^0 
could be taken. We present the results as the sum of two contributions. For the 
first we will add together the contributions from sections 3 and 4 with a factor 
—7r^/3 for both. In the second contribution, we compensate by using the result of 
section 4 but with a factor —tt^/G. 

5.1. Contributions from Figs.7, 8. We take the former contribution first, com¬ 
bining all the terms arising from the graphs in Figs.7 and 8. Since we have complete 
contributions (not separating Fig.8 from Fig.7) the ghost graphs disappear in the 
limit 0 —>■ 0, because of the factors [1 — (P^/P^)] in Fig.5. Thus the result comes 
from graphs (a), (b), (c) in Fig.7 and (a), (c) in Fig.8. It is straightforward to write 
these out using the rules in Figs.3 and 5. 
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Fig. 13 Graphs with n-1 denominators contributing to the foctorB in 
equation (11). They are derived from Fig.7 using Fig.5. 


We first note that, with 0 = 0, the only external momentum dependence in (11) 
comes from the second 5 function:, which we may write 


(29) 


HPo + Qo+ro) = ^ 



L 771 



where the (m = 1, 2, 3...) are the momenta in the external fields. The effect of 
this is that the external field A“(K, ko) is everywhere replaced by the time Fourier 
transform A°;{K,t). From now on, that external fields in (7), (9) and (31), (32), 
(33) below are to be understood as meaning this. 

Then the contributions have the form 


(30) 

-(gV8)(27r)-^+V“'’"r''’'^' [ dt [ dPdQdRdP'dQ'dR'S{P+q+Tti)S{P'+Q'+Il')F, 


where the terms in F come from the graphs mentioned above. (In (30), the integrals 
and delta functions are understood to be (d—l)-dimensional.) Expression (30) is for 
the effective action; for an effective Hamiltonian, the time integral and the minus 
sign should be omitted. 

Define 


j“/(p,p') = r,,/(P)(5““'d(p - p') 
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(31) + J dP"Ti»i(P)yl““"(P-P")G“"“'(P",P')(zP'-) 

and 

(32) L““'(P,P')= ydP"T,i»(P)7l“>(P-P")G“"“'(P",P') 

where the transverse projection operator Tiiii(P) = 5i"i — {PiPin)/P"^) comes from 
Fig.5 when 0 = 0. With this notation Fig.7(a),(b),(c) together with Fig. 8 (a),(c) 
give a contribution to F in (30) 

(l/3)[-J-ih-P', -P)G“'(R, R')4''(Q, Q') 

+J“'“(-P', -P)G“'(R, R')if (Q, 

(33) +(7R,)Lr(-P',-P)G='^(-R',-R)jf (Q,Q')]. 

5.2. Contributions from Figs.15, 17. In this case, unlike the previous one, even 
when 0 = 0 , we cannot neglect ghost graphs, because only incomplete ghost graphs 
appear in Fig.17 (the 1 parts from Fig.5 without the —P^/P^ parts); so we must 
take all the graphs. As explained above, they are taken with a coefficient 1/6. 
Thus the contributions have the form 

(34) (l/ 6 )[(ci)+(d 7 )'+(ei)+(da)+(cii)+(eii)] 
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Q = kl toplM 
[] = 1 in bottom ine 


Fig,16 Graphs with r derominators obtaired iionn Fig.8. Theq-irie(d[awninltie<»i^Bs 
a IMn rectangla) iias& In (a) ard (c), and 1 In {d) and (a). Tbesaoond ldan% In 
Fig 6 haa tiaan usad (Mnca tiMO conbrlbutlQnsItiini (e)and (a])i. 



Fig.'l 6 Ttiroo graphs obtained from Fig.'l S (d) by moinsrhan oonservatior. 


giving 

(35) (l/3)[(ci)+(ei) + (cii) + (eii)]-(l/ 6 )[(cii)+(eii)]. 

For the first square bracket in (35), we use the identities in Figs.5 ,6 to get back to 
graphs like (c) and (e) in Fig. 8 , but with the g-line replaced by Sij and 1 respectively. 
These are shown in Fig.19, and have a coefiicient (1/3). This must be doubled 
because of graphs with the p- and ( 7 -lines intercanged. The result is (with the 
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Fig.17 Graphs with n-1 denominators obtained from Fig.8. ITw q-line 
(drawn in the center as a thin rectangle) has just a or (in (d)) 
a 1. The notation of Fig.5 is used. 


notation of (32)) 

(36) (1/3)[2L“'“(-P', -P)][5“',5(Q - Q')][*i?»G^'^(-R', -R) - ii?'G=='(R, R')]- 

There is also the contribution from the second square bracket in (35), with 
coefficient (—1/6), but there is also a minus sign with (cii) and (eii) in Fig.15. This 
also must be doubled, giving 

(37) (l/3)[iP,G““'(-P',-P)][<5“'(5(Q-Q')][ii?*G^'‘=(-R',-R)-ii?'G“'(R,R')]. 

The expressions in (33), (36) and (37), inserted for F in (30), constitute the 
generalization of the Christ-Lee terms when the external field is not transverse. 

6. The transverse special case 

In order to make contact with 0, [6], we show how our result can be re-written 
if the external gluon field is transverse. In this case, many simplifications occur. 
The covariant derivative D obeys d.D = D.d and so 

(38) G“'(P,P') = G“'“(-P',-P), 

that is the arrows in Figs.3 and 6 and subsequent figures are not necessary. 

We examine the contributions to F in equation (30) in this special case. We note 
that in Fig.17(d), instead of using the Jacobi identity as we did to get Fig.18, we 
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+ 0 



(dl) 






i8> 


F1g.1 fi Graphs oblalrad fram Rg.17 (d) by da jbla usa of Jacohl Idarlity 



Fig.19 cartitjutiig to tba fmction F in equation (30), which are derived fiom 

Rgs 15, 16 and 17. They have caeffldanlB V3. The ramalnlrg corblbutkins 
fram Fig. 15 are the tiwD graphs (dl) and (ell), which have coentdart -1 /6. 


0 








Rg.20 Graphs obtaired from Fig.17 (d) by using momentum CGmervahan. These oe 
relevert to the trenaveree spedel cnee, eo the anroee ere orriilted ftoiri the left 
hand Unas. 


might alternatively use momentum conservation, which yields Fig.20. We will take 
the average of these two; that is Fig.18 with coefficient 1/6 together with Fig.20 
with coefficient 1/6. Each of these contributions has to be doubled because of the 
graphs with p, q interchanged. 

Then we take 


(l/3)[Fig.l8 + Fig.l5{(ci) + (ei)}]+(l/3)[Fig.20 + Fig.l5{(ci) + (ei)+(cii)+(eii)}] 



QCD EFFECTIVE ACTION 


21 


(39) = (l/3)Fig.l9 + (l/3)Fig.21, 

where we have used the first and second identities in Fig.6. (We neglect terms like 
(28) which are zero in dimensional regularization.) 

In [B], equation (4.5.13), it is shown that the contributions from Fig.7(a),(b),(c) 
and Fig.8(a),(c) together with (39) are equal to Vi + V 2 of [S]. Thus we have 
confirmed that, in the transverse special case, our method agrees with the known 
results. 



Fig Graptis corvbijlhg b tlie flncbor F In equatton (30) derived fbri Flge.15,16, 
17 and 20. Thsss are relevari b tin bersveiee cbbs, so ttie anows are oirhted 
on the left. 


7. SUMMARY AND DISCUSSION 

We have proved that, in spite of energy divergences in individual 2-loop graphs, 
the effective action in perturbative Coulomb gauge QCD is well defined. The sum of 
the relevant terms (in expressions (33), (36) and (37)), which is energy-independent, 
is calculated, being a generalization of the function Vi -I- V 2 derived in a different 
way in [ 5 . These terms in the effective action would be a necessary contribution 
to the BRST identities. As a check, we have shown that in the special case of 
transverse external fields, our result reduces to the known one. 

What are the implications of the Christ-Lee function Vi + V 2 and its gener¬ 
alization in this paper? There are two quite different approaches, the Christ-Lee 
approach and the Doust approach . The Former carefully studies the operator 
ordering in the Coulomb Hamitonian, and derives the Hamiltonian 

(40) Hw + Vi+ 14 

where the suffix W denotes Weyl ordering, which is necessary to avoid double 
counting. This is because of the zero vacuum expection value of 

(41) {AU'P)E^,{P)}w = (l/2){i“(P)A^(P) + A^(P)i“(P)}, 

which implies that Feynman graphs like Fig.l, with pairs of operators contracted 
in the Coulomb potential operator, are omitted. 

The Doust approach, followed in this paper, effectively defines the Coulomb 
gauge as the limit of the flow gauge given by equation (4), and there is no operator 
ordering problem in the flow gauge. 

Presumably the Christ-Lee Hamiltonian (40) might be used in Hamiltonian lat¬ 
tice QCD, but the Doust approach would not be appropriate. On the other hand. 
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the effective action for non-transverse fields cannot be derived in the Christ-Lee ap¬ 
proach; and this general effective action is required to implement BRST identities. 

Appendix - Definition of effective action 

The effective action is normally defined by first introducing a source current J: 

(42) J d<^xJt{x)A<t{x) 

(for simplicity, we omit sources for other fields). Then the external classical field 
A1(x) is introduced by a Legendre transformation from J to A. If the transversality 
constraint is imposed by a delta function: 

(43) 5{A<t] 

there appears to be a problem, because then only transverse Jf contribute in (42). 
To avoid this, we may impose transversality by using a gauge-fixing term like 

(44) ^ I d^^xid.Atix)}^ 

and taking the limit 0 —>■ 0 only after the effective action has been defined. Actually, 
for the purpose of the present paper, we should use the gauge fixing term in equation 
(4) in place of (44). 
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